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ABSTRACT 


A method for obtaining the non-geostrophic component of the 
observed wind is outlined. Two techniques are used to obtain the | 
divergence and relative vorticity of the non-geostrophic component. 
The effects of several influences on the ageostrophic divergence and 
relative vorticity are discussed, including condensational heating, 
gradient winds, and terrain effects. 

Two case studies of the 700-400 mb layer over a portion 
of the Western Cordillera are presented. The results of these studies 
indicate that horizontal thermal advection was the dominant influence 
controlling the ageostrophic wind, with condensational heating and 
gradient wind influences having only minor effects. It is concluded 
that prevailing atmospheric stability and the methods of analysis 
employed did not permit sufficient resolution of the wind field to 
delineate terrain influences in the cases studied. Some suggestions 
for further work are made within the scope of the methods discussed 


which may help to determine terrain effects. 
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CHAPTER I 


INTRODUCTORY REMARKS AND THE CALCULATION 
OF DIVERGENCE AND VORTICITY 


fer Lutrodicl ron 

Approximations are an integral part of meteorology. Without 
them, most problems would be so complex as to deny solution. However, 
the pressure on the meteorological community to produce better forecasts 
covering longer periods of time has forced meteorologists to re-examine 
the ramifications of the approximations used in dealing with atmospheric 
phenomena. Interactions between the Earth's surface and the atmosphere 
have emerged as one of the major areas where rough approximations are 
no longer sufficient. In particular the Western Cordillera is a factor 
which must be given due consideration: it affects atmospheric motions 
on all scales -- from the valley wind to the planetary wave. The advent 
of the digital computer and the profusion of numerical circulation 
models which followed, have done much to clarify the mountain-atmosphere 
interaction. This clarification (for synoptic scale features) has come 
about principally through the use of the @-equation in conjunction with 
the quasi-geostrophic equations (cf. Haltiner, 1971). As Sawyer (1959) 
points out, the quasi-geostrophic system does not permit the evaluation 
of ageostrophic effects since this type of motion is ignored by the 
equations. While this does not appear to be a serious defect for broad 
(“1000 km.) and low (maximum height ~ 300 m) mountains, ageostrophic 
effects arising from a barrier having the dimensions and complexity of 
the Western Cordillera may be significant. 

On the synoptic scale, terrain-induced vertical motions have 
an important influence on developments in the vicinity of mountains. 
Since these induced vertical velocities must eventually be damped out at 
some level above the terrain, it follows that they must create a 
localized field of horizontal divergence and convergence. Namias and 
Clapp (1946) discovered this field while using monthly inean maps of the 
10,000 foot level ( «700 mb). They remarked on its seasonal consistency 
and attributed it to topographical effects. To the extent that the 


ageostrophic wind is neglected in the quasi-geostrophic numerical 
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formulations, it follows that the divergence and relative vorticity 
of this component will also be neglected. It is the intended purpose 
of this study to try and determine if the ageostrophic wind component 
over the Western Cordillera shows any evidence on a synoptic scale of 
spatial organization that can be attributed to terrain influence in 
terms of divergence and the vertical component of the relative vorticity. 
The data-gathering network around and in the Cordillera is 
not sufficiently dense to permit an investigation at any level below 
that of the highest level of the terrain i.e. »3000 m in the region 
studied. The 700 mb level was chosen as the lowest acceptable level to 
be examined. Induced vertical motions will be the principal agent by 
which terrain effects are communicated to the upper levels of the 
atmosphere. At the 400 mb level these velocities will generally be 
about one order of magnitude less than their initial value (Berkofsky, 
1964), hence the analysis was terminated at this level. 
The methods employed will be simple in view of the lack of 
an adequate theory dealing with ageostrophic winds. The ageostrophic 
components will be found by subtracting an objectively-analyzed 
geostrophic wind from the observed wind vectors. The calculation of 
the divergence and vertical component of relative vorticity is 
discussed in section 1.2. Later discussion (section 2.5) will show 
the limited accuracy of these methods, with the result that the presence, 
and not the magnitude, of ageostrophic effects will be of prime concern. 
Bt should tbe noted at this point) that the @x;y5p;t) 
coordinate system will be used throughout this thesis. The coordinates 
x and y will refer to the finite difference grid in the normal Cartesian 
sense. In this respect uzdx/dt, and vedy/dt. All motions will be with 
reference to a constant pressure surface, and subscript "A' will denote 


ageostrophic. 


1.2 The Methods of Divergence and Vorticity Calculation 


To obtain consistent and reasonable estimates of atmospheric 
vorticity and divergence is a continuing problem in synoptic 
meteorology. Numerous techniques have been proposed and used for 


computing these quantities. The method chosen depends on many factors, 
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some of which are: 

-- The purpose for which the quantities are being obtained; 

-- The required accuracy; 

-- The quality and amount of initial data required or 

available; 

-- The ease of computation; 

-- The atmospheric processes to be taken into consideration. 
Deciding on the method to be used usually involves a compromise of one 
or more of these factors. Two methods appropriate for the ageostrophic 
winds were chosen for use in this study. The first method will be 


called the Fundamental Method in which 


Glv2 ed) 


where D, and Qe are respectively the horizontal divergence and vertical 


F 
component of relative vorticity in the (x,y,p,t) coordinate system. 
The principal advantage in using (1.2.1) lies in the ease with which 
these equations can be evaluated. However, the fundamental method has 
an inherent disadvantage which must be recognized when it is used. 
This drawback is its sensitivity to errors and Gaussian noise in the 
initial observations, and in subsequent analyses of both the pressure- 
surface heights and the wind vector. Observational errors can be 
partially removed by careful checking, but Gaussian noise is difficult 
to remove because of the rather arbitrary nature of the choice as to 
what constitutes noise. While methods of polynomial fitting and 
harmonic analysis have been devised for minimizing the above mentioned 
errors, their usage is again limited by what is chosen as the level 

of significance. 

As a compromise, the decision was made to use the time- 
averaged ageostrophic components as the basis for calculating the 
divergence and vorticity by the use of (1.2.1) under the assumption 
that the actual uncertainties in the values for specific times would 
tend to cancel in the averaging procedure, even though the theoretical 


uncertainty is increased (section 2.5). The averaging may also remove 
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or minimize unwanted ageostrophic effects accompanying minor 
disturbances moving through the region studied. The shifting of the 
flow patterns during the interval over which the averages were taken 
will mean that the calculated results will apply to the mean flow of 
the period. 

The second method of estimating the vorticity and divergence 
of the ageostrophic deviation is that proposed by Bellamy (1949). If 
it can be assumed that variations in the wind field are linear between 
any two observing stations forming the vertices of a triangle, then 
it can be shown (Appendix I) that the divergence and relative vorticity 


of the flow contained within the triangle are given by: 
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in which subscript i represents the oo vertex of the triangle, Va is 
the ageostrophic wind speed, p; is the azimuth of the opposing side 
and «, is the direction of un (both measured clockwise from the 
positive y axis of the grid), and h,is the length of the perpendicular 
to the opposing side through the i vertex, 

As employed in this study, the Bellamy method requires 
Several considerations. First, since it was used on the ageostrophic 
wind field, it will be sensitive to the accuracy of the observed 
winds and pressure-surface heights. Second, the triangles employed 
should be as close to being equilateral as possible. This is necessary 
to ensure that the results represent two-dimensional phenomena. To 
satisfy this condition, no triangle should have any interior angles 
greater than 90°. 

The third factor to consider is the point in space to 
which the Bellamy estimates apply. If the triangles were perfectly 
equilateral, the obvious choice would be the point at which the 
bisectors of the interior angles intersect. The triangles used in 
this study were not truly equilateral, with the result that the 


"center of gravity" of the triangle was chosen as the point to which 
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the results would apply. Eddy (1973) indicates that this will introduce 
a geometrical bias into the results. 

The assumption of linear variations in the ageostrophic 
wind field between stations appears reasonable with the scale of 


analysis used and the size of the variations involved. 
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CHAPTER II 
POSSIBLE SOURCES OF THE AGEOSTROPHIC WIND DEVIATIONS 


2 ie inerodwetion 


To obtain the ageostrophic deviation, the wind field was 
=> 


decomposed into two components -- one in geostrophic balance ( V 


3 


and a non-geostrophic component ( Wa ). Hence, the observed wind 


—_—> 
( V) can be represented mathematically as: 
V = Mai Nie (Fe 41) 


thus: from 24.) )iwethave: 


Uw, = U- Ug 
C2oh,2) 
Vp = V- Vy 
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v ge 9(9%) 

f\dIx/p 
in which x is the geopotential height of the reference surface p, 
the acceleration due to gravity is g, and f represents the Coriolis 
parameter. 

To presuppose that the ageostrophic deviations calculated 
are generated solely by terrain influences would be folly since 
ageostrophic deviations may arise from many factors. However, to 
try and decompose them into components due to each, while theoretically 
possible in some cases, would be very difficult for practical purposes. 
In lieu of this, the rest of this chapter will be devoted to discussing 
the effects that the influences listed below would have on the values 


of divergence and vorticity of the ageostrophic wind, to aid in 
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interpreting the results: 
(i) Atmospheric compressibility, 


"oradient wind") 


(ii) The presence of curved flow (i.e. the 
(iii) Thermal influences due to condensation processes, 
(iv) Propagation of measurement errors, 


(v) Terrain influences. 


2.2 The Assumption of Incompressibility 

Batchelor (1967) employed the methods of scale analysis to 
the study of the conditions under which the assumption of incompressib- 
ility was justified when considering the velocity distribution in a 
moving fluid. He assumed that the fluid flow can be characterized by 
a scale length L such that variations in the fluid velocity Vv, are 
small over distances that are small compared with L. Furthermore, 
the variations of VY in both time and space are assumed to have the 
magnitude U . The order of magnitude of the spatial derivatives 
of components of Vv will then be given by U/t.. For the purposes 
of this study, the length scale was chosen to be half the width of 
the main current at 500 mb which was found to be about seven degrees 
of latitude, or roughly 8x10° meters at 45°N. The scale velocity 
thus becomes the central value of this current which has a typical 
value of about 40 m eeeme. 

Let us assume an atmosphere having constant composition 
and choose as the parameters of state the density P and specific 


entropy S. Then the rate of change of pressure experienced by a 


material element can be expressed as: 
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Rearranging (2.2.1) gives: 
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and it is to be shown that 
ive Vicou7 ls 


Hence the assumption of incompressibility will be justified if: 
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At this point the mathematics becomes somewhat involved. 
The total derivatives in (2.2.3) are expanded (using the isentropic 
assumption to expand dp/dt) and the condition is imposed that all 
terms must be much less than Ai The reader is referred to 
Bachelor's text for the details, but the resulting inequalities are 
given by: 
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in which: 
es ; F ° 
WY = A measure of the dominant frequency of oscillation 
in the flow field, chosen to be every three days or 


about 4 x 107° ee 


~® = The isobaric coefficient of thermal expansion 


= 3.92 x 107° °c’ @ -17.8°C and 1013 mb 


(2) 
~O 
“" 


Specific heat at constant pressure of dry air 


=) 40. 240 cal. an ce 


= 


Vy the kinematic viscosity 


= 0,132 cm’ sec’ @ 0°C and 1013 mb 


@ = A measure of the temperature differences in the fluid, 


chosen to be 50°C 


Ww = The thermal diffusivity of dry air 


= 0.184 om” sec | @ O°C 


C = The local adiabatic speed of sound 
320 m sec @ 5 km 


ree 


The scale of vertical motions, chosen to be the 
height of a homogenous ( P = constant) atmosphere 


or is approximately 8 x 10° m 


Wale Magnitude of vertical velocities, chosen to be 


approximately 0.1m et: 


Substitution into the inequalities (i)-(v) shows that 
all left-hand sides are at least two orders of magnitude less than 
unity. This was considered sufficient to warrant neglecting atmospheric 
compressibility, and hence its contribution to the divergence or 


vorticity of the ageostrophic wind. 
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2.3 Results Arising from the Gradient Wind 


The flow patterns of the atmosphere, in the cases chosen 
for study, show evidence that in some areas the consideration of a 
gradient wind would be appropriate. Unfortunately, the necessity 
of calculating the curvature of the flow makes the gradient wind a 
difficult quantity to calculate with any degree of accuracy. 


The definition of a streamline may be expressed as: 


CAE &) 


EIR 


“4 = 


where y is a function of x only. If we assume that the observed wind 


speed (V) is in fact the gradient wind speed G, then: 


We aig ely (2.3.2) 


Let us further assume that the flow can be represented by an 
infinitely broad current flowing from west to east in which the 


equation of a streamline is given by: 
> a Cees (23353) 


where a,A and “AQ are constants and /\ can be any desired contour 
spacing. A Value of a=l and } =20 was found to yield a contour shape 
shwon in Figure 1 which is roughly similar to those of interest in 


the case studies if x and y are in units of grid-interval. 
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Petterssen also gives the following expressions for the 


VOreLeLtLy | Q ) and divergence ( J) ) of the wind along a streamline: 
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(2.3.4) 


where p is the angle which the wind makes with the x axis, Ky is the 


curvature of the streamline, and s and n meaSure distance along and 


perpendicular to the left of the streamline if moving in the direction of the 
flow. For the purposes of this analysis, the streamlines, trajectories, 


and contours will be assumed to be coincident and given by (2.3.3). 
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Fig. 1 The hypothetical field of flow for estimating the gradient 
wind contribution to the ageostrophic divergence and relative vorticity. 
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Substitution into (2.3.4) gives: 


Peres ee 
2 98 is 
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where the substitution K =1/R has been made, R being the radius of 
curvature of the contour. Subscript G denotes "of the gradient wind". 

The final assumption to be made is that the gradient wind 
speed (G) involved in the equations (2.3.5) is the maximum possible 
at any point along the contour. Making this assumption, and setting 
G to some constant value such that it will not be a function of 
position, will lead to calculated results which are the maximum 
possible at any given point along the contour for the chosen gradient 
wind speed. 

Using these assumptions, equations (2.3.5) were evaluated 
(see Appendix I) for the flow pattern of Figure 1, assuming a value 
for G of 40 m yeep? the results are shown in Figure 2. 
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Fig. 2 The divergence (dotted line) and vorticity (solid line) of 
the gradient wind for the flow field of Fig. 1 using G=40 m sec-l, 
Ordinate units are sec-!, abscissa units are taken to be grid intervals 
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An examination of these results indicates that the relative 
vorticity of the gradient wind may be of the order of to Boca. for 
strong flow exhibiting moderate curvature. The divergence arising from 
this wind does not appear to be significant, having a maximum value of 

-6 = hi -1 ; 
order 10 sec ., For values of G other than 40 m sec , for instance 
G', the relationships: 


De (8)Ps 


re 
Ry -($5) on 


can be employed using the results in Figure 2. 
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2.4 Thermally Induced Vorticity and Divergence 


Danard (1964) in a case study of a rapidly developing cyclone 
found that the release of latent heat due to condensation in the free 
atmosphere could result in considerable changes in the large-scale 
vertical motions. Vertical velocities in areas of active condensation 
were found to vary by a factor of two or three from those calculated 
when release of latent heat was neglected. In a previous but similar 
study by Aubert (1957), the same conclusiins were drawn, and it was 
shown that these effects are rapidly damped as one moves away from the 
region in which condensation is occurring. 

These two studies indicate that a consideration of thermal 
effects can considerably alter the derived patterns and magnitudes 


of vertical velocity. Thus, since: 


pe eee ids (22719 


where D_ represents the isobaric divergence on a constant pressure 

P 
surface and a) = dp/dt, itis evident that an inclusion of thermal 
effects may alter the magnitudes and distributions of divergence 


and vorticity along a constant-pressure surface. 
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Let us assume that the wind speed at a point in the free 


atmosphere, with no condensation occurring, is given by: 
Wee Vy = Wy k (2.4.2) 


Let us further assume that if condensation occurred at this 
point, then the vertical speed would be doubled, while the total wind 


——— 
( Ne ) remained the same such that: 


= 


l= Vet aw, k (7s) 


—— 


where Ne represents the wind components along the pressure surface 
with condensation. Solving (2.4.2) and (2.4.3) for the change in 


horizontal wind speed gives: 


~N 
ae W, \ 
where 
es AEE eek, 
= NE | 
Da ue ve 


Thus, in order for the horizontal wind to change by only 
one meter per second, the initial vertical wind speed would have to be 
at least one meter per second. Vertical velocities of this magnitude 
generally occur only in the most intense subsynoptic systems. A more 
likely vertical velocity of say 10 cm Beye ts would require an adjustment 
or ronly20.1 m cae in the horizontal components. Resolution of the 
wind vector to this degree with the data available is not possible. 

In view of this, quantitative estimates of the effects of 
condensational heating were considered impractical with the methods 
of calculation employed. However, qualitative judgments were made if 


such were considered necessary. 
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2.5 Propagation of Errors Leading to Spurious and/or Uncertain Results 
In as much as this study is basically an analysis of 
atmospheric measurements, it was deemed necessary to examine how the 
errors in these measurements would affect the calculated values of 
divergence and vorticity. The method of error analysis chosen was 
that’ outlined "by ‘Baird "(1962)4 
As part of the data check, the heights of the 700, 600, 
500, and 400 mb surfaces computed using the hydrostatic assumption 
were compared to the values reported on the ascent recordings. If 
the observed heights are accepted as the correct values, then the 


probable error (P.E.) in the calculated values would be given by: 


‘P ae Gan (2.5.1) 


where omy represents the variance of the difference in the 
two height values. 

However, the observed height values are not the "true" 
values, but incorporate errors from several sources, such as instrument 
and observer errors. An estimate of the observational errors was 
obtained from a study done by Rapp (1952). The probable error 
calculated by the method given above, and Rapp's estimates of the 


error in the observed heights are given in Table l. 


Pressure Level (mb) PE. yf Spm) P.E. oy (spm) 


700 Lae = 4.7 
600 pel ee 
500 2.4 6.8 
400 es, ‘ala 


Table 1. A comparison of the probable error in the 
calculated heights accepting the observed value as correct (P.E. 
and the probable error in the observed heights (PE e) < 
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Since the difference between the calculated and observed 
heights lies well within observational error, it was decided that the 
observational errors given by Rapp could be used as the uncertainty 
in the calculated heights ( S¢ ) of the various pressure surfaces. 

The uncertainty in the heights of the pressure surfaces 
will lead to an uncertainty in the calculated geostrophic winds as 
defined by (2.1.3), which in turn will lead to further uncertainty 
in the ageostrophic components defined by (2.1.2). Note that there 
will also be a contribution to the error in the ageostrophic components 
due to uncertainties in the observed wind speeds OSES Rapp's 
values for the probable error in the observed wind speeds are given 


rin MAI: 2. 


Height (km) PLE. (m cea} PLE. (m ae 
0.5 On21 0.35 
EXO 0.45 0.45 
LD 0.38 0.27 
220 0242 0.18 
ots OFS) 0753 
3.0 Of43 0.41 
Shas) 0.24 O19 
4.0 Dose 0.60 
4.5 0.46 0765 
Se0 0.83 0.99 
Ae, 0.90 0.69 
6.0 0.70 1.00 
655 0.49 50 
120 0.40 0.54 
yee O23 1.80 
8.0 0.49 0.84 
8.5 1.99 2500 


Table 2. The probable error in the observed wind components 
Calter nappy. Looe). 
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Rapp attributes the inconsistency in the P.E. of the 
wind speed between consecutive levels to the sampling technique 
employed. Since the errors at all levels are of the same magnitude, 
an arithmetic mean was used for the values from 0.5 to 7.0 km 
inclusive, yielding values of SUW=0.48 m sec” and SU =0.60 m secre 

It should be noted that Rapp's error estimates are 
probably very conservative since the observers in the experiment 
and the conditions under which the raw data were collected were of 
a higher quality than those involved in the routine program of upper 
air sounding. Carrying through the error analysis as given in 
Appendix III gave the following expressions for the uncertainties in 
the divergence and vorticity as calculated by the fundamental 


(subscript F) and Bellamy (subscript B) methods: 
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5 Qoz *a/ —Ua__)( Surgst (0.5.4) 


An examination of the ratios SA and Ory fay Oe 


shows that the results of the Bellamy method will generally be less 
uncertain than those of the fundamental method. Thus, if the two 
methods yield similar results, then some semblance of the actual 
physical situation can be attributed to the calculated values 
inspite of the large potential for error shown by the sample results 
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Level (mb) Sp.» 5Q, (see) 
700 ates On. 
600 37 
500 42 
400 67 


Table 3. The uncertainties in the ageostrophic ass 
and vorticity given by (2.5.2) using H=185.3 km, §v=0.60 m sec +, 
fel.1xl074 sec-l at 50°N, and Table 1. 


2.6 Divergence and Vorticity Induced in the Upper Flow by Mountainous 
Terrain 
If an air stream impinges on a mountain barrier, the terrain- 
induced vertical motions will have an effect on the air stream. 
Berkofsky (1964) studied three methods of computing the fall-off with 
height of the terrain-induced vertical motions. He suggested that 


equation (2.6.3) is a viable alternative to more complicated methods. 


ads (2) (4x) (2.6.3) 


Where (Jp, is the vertical velocity at the initial level 
po and: 
pec il eel 
peep 
is the stability factor with: 
n= R = 0.288 
Cp 
Studies by Charette (1971) and Hopkinson (1972) suggest 
that terrain induced vertical motions over the Western Cordillera can 
be of the order of +10-20 cm Rao (or roughly +LO-20 7) eee) at 
the 850-mb level. Using these figures and the radiosonde data for 


the case studies together with equations (2.4.1) and (2.6.3) indicated 
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that terrain induced divergences could be of the order of 

ty 0=30xK10 F coche at the 625-mb level. While this is a very crude 
estimate, it does serve to indicate that in regions where extensive, 
and abrupt changes in topography occur, such as the windward and the 
lee slopes of the Western Cordillera, the resulting vertical motions 
may generate divergence at levels well above the terrain that is 
synoptically significant if the motions are sustained over a large 
enough area, 

No suitable method could be found for estimating terrain 
induced vorticities in a simple fashion. Clearly, the largest terrain 
effects must occur at the point at which the atmosphere and the surface 
are in contact (i.e. the boundary layer). If one assumes that the 
terrain itself is a pressure surface such that variations in the 
height Z of the pressure surface are given by variations in the 
height of the terrain, then the expression for the relative vorticity 
Q of this hypothetical surface assuming the Coriolis parameter f does 
not vary would be given by: 

2 

Q= 9, Vp h (2.6.4) 

The values of Q so produced are of the order of Om cece 
assuming a grid interval of 185.3 km in the area of interest (see 
Section 3.1). While the processes of atmospheric dissipation rapidly 
reduce this value, some portion must be transferred to higher levels 
of the atmosphere by vertical motions. However, the nature of the 
decay is not certain; hence the portion of the calculated ageostrophic 
vorticity arising from terrain effects must also be uncertain. In 
view of these problems, the calculated relative vorticities will be 
compared with the terrain Laplacian for agreement in sign, with due 


allowance given to the other factors discussed in this chapter. 
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Summary 
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Having shown in Section 2.2 that the assumption 

of incompressibility is justified for the scales of 
motion involved in this study, divergence and vorticity 
arising from compressibility effects were considered 
undetectable, and therefore neglected. 

The contribution to the ageostrophic wind from 
condensational heating is felt to be undetectable 

with the quality of the data available. However, the 
effects on the magnitude of the calculated divergence 
may be significant and will be acknowledged in a 
qualitative manner. 

The flow patterns indicate that the cyclostrophic 
component of the gradient wind may contribute 
significantly to the ageostrophic wind. In areas where 
this occurs, allowance was made in the calculated 
divergence and vorticity by using the results in 
Pieure 2 and (223.60). 

Propagation of observational errors through the 
calculations leads to large uncertainties in the 
divergence and vorticity. A rough estimate shows that 
the values calculated by the Bellamy method will generally 
be less uncertain than those of the fundamental method. 
In interpreting the calculated fields-of divergence 

and vorticity, more emphasis will be placed on those 
regions where the results of the two methods are 
comparable. 

Terrain induced vertical velocities may make a 
significant contribution to the field of horizontal 
velocity divergence. The magnitude of this contribution 
will depend on the speed and direction of the low level 
flow (i.e. below 700 mb) since this will affect the 
initial vertical velocities. However, a rough estimate 
indicates that terrain induced divergence of the order 


of 10-30x107° Eaeae may exist. The terrain contributions 
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to the ageostrophic vorticity patterns at and above 
700 mb are questionable, hence the patterns will only 
be compared with the terrain Laplacian for qualitative 
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CHAPTER IIL 


THE DATA AND THE ANALYSIS PROCEDURES 


3.1 The Finite-Differencing and Objective-Analysis Grid 


The grid employed was designed to be used in conjunction with 
a polar-stereographic map projection having a standard latitude of 60°N. 
The grid interval H used is 185.3 km. The reference meridian chosen 
is 80°W since it is roughly parallel to the Alberta portion of the 
Rocky Mountains and facilitates the interpretation of the results. An 
initial grid of 14x10 grid intervals was used giving 165 grid points. 
As shown in Figure 3 the geographical area covered includes Washington, 
Oregon, Idaho, and most of British Columbia, Alberta, Wyoming, Montana, 
and Nevada. The radiosonde stations used for the raw data are shown 
im Ergure 3. 

The equations for finding the grid coordinates of the data 
points from their latitude and longitude are derived in Appendix IV. 


The pertinent two relations are: 


yg = 15 ty ( tan Yr cos xp—tan ¥s Cos 3) (351.19 


2 A 


X= eal ( tan Sinx gs — tan *¥p sin <p) (4.192) 


H pl 2 


in which subscripts p and s represent "principal-grid-point" and 


"station" respectively, WY is the co-latitude, Sy ae the angle 


? 


between the reference meridian (80°W) and the meridian through the 
point, J=R,(1+cos YW. ) where Re is the radius of the earth and W is 


the standard co-latitude. 
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Fig. 3 The finite-difference and objective-analysis grid. Stations 
shown are the data points. The thin dotted line is the approximate 
position of the 1500 m ASL terrain contour, Grid interval is 185.3 km. 
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Fig. 4 The geometry for transforming the observed wind direction 
to yield wind components which are referred to grid north. 


The observed wind direction must be transformed in order 


to obtain component velocities which are consistent with the grid 


employed. Using the geometry shown in Figure 4, the transformed wind 


components are given by equations (3.1.3) 


Cee Ni sin ¥ 


(Soi.2) 
ae gC 0 at 
where %= Be SG if A-% 37560 

¥=0 if B-%,=360 

$= A- Xs ie Os A eae eo 


and 6 is the wind direction measured clockwise from true north. 
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3.2 Synoptic Description of Case Studies 

Two cases were chosen for this study. The first case 
covered the period from 21 December 1200 GMT to 24 December 0000 GMT 1971, 
inclusive. The upper-level flow during this period was dominated by 
a strong low in the Gulf of Alaska which moved rapidly out of the Gulf 
at the end of the period to a position about 300 nautical miles west 
of the Oregon coast. A weak quasi-stationary ridge from Central 
Alberta to the Dakotas persisted through the period and maintained 
a strong to moderate flow from the southwest over southern British 
Columbia and Washington State. Dewpoint depressions over the area 
at 700 mb were in the 0-5°C range for most of the period. 

The surface charts showed a double-centered low pressure 
system, with one center near Queen Charlotte Island, and the other 
about 100 nautical miles southwest of Vancouver Island. This system 
was bringing maritime air into southern British Columbia and Alberta 
as well as the Pacific Northwest states. An Arctic high was over 
Alaska and pushing continental Arctic air into Alberta from the northeast. 
During the period the Arctic high broke down, but an unstructured 
region of high pressure was set up over Alberta and Saskatchewan. The 
low centers combined and shifted slightly south to join with a low 
which moved up the United States coast. The baroclinic zone shifted 
from southern British Columbia and Alberta into the Washington, Idaho, 
and Montana area. 

The second case studied in detail covered the period 
16 January 1200 GMT to 18 January 0000 GMT, 1972 inclusive. The 
Synoptic features were similar to the previous case. At the beginning 
of the period the upper flow was characterized by a strong low over 
the Gulf of Alaska with a trough extending into northeastern 
Saskatchewan. During the period the upper low weakened into a broad 
trough on a NNE-SSW line through northwestern British Columbia. A 
weak ridge from Edmonton to Great Falls, Montana, persisted through 
most of the period. The major current shifted south from central 
British Columbia and Alberta to a position near the United States 
border, maintaining a rough east-west orientation throughout the 
period. Dewpoint depressions at 700 mb were in the Q-5°C range over 


most of the region for most of the period. 
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Surface features were marked by a strong low in the Gulf 
of Alaska bringing maritime air into southern British Columbia and 
Alberta. During this period an Arctic high moved into northern 
Alaska and stalled there, pushing continental Arctic air into Alberta, 
and driving the baroclinic zone south through Central British Columbia 
and Alberta. By the end of the period the Gulf low had weakened 
considerably, but exhibited only a minor southward shift, while the 
baroclinic zone extended through southern British Columbia and Alberta 

In both cases studied, there was considerable cyclonic 


activity along the baroclinic zone although no major developments 


ensued. The regions over which continuous precipitation occurred during 


the periods was analyzed and an average precipitation area produced 


for each case as shown in Figures 5 and 6. 
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Fig. 5 The average areas where precipitation occurred during the 
period 21 December 1200 GMT to 24 December 0000 GMT, 1971, over that 
portion of the grid where the ageostrophic winds were analyzed. 
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Fig. 6 Same as Fig. 5 except for the period 16 January 1200 GMT 
to 18 January 0000 GMT, 1972. 


3.3 The Aerological and Terrain Data 


Since analyses were desired at levels not reported in the 
monthly summaries of radiosonde ascents, it was necessary to obtain 
the records of complete upper-air soundings. Data for Canadian 
radiosonde stations were available from the Atmospheric Environment 
service as Adiabatic Charts and Winds Aloft Computation forms. Data 
for United States stations came from the U.S. Department of Commerce 
Environmental Data Service as WBAN 31 Adiabatic Charts, WBAN 20 
Winds Aloft Computation Sheets and RAOB Computer listings. All data 
were obtained on microfilm. 

Pressure, temperature, relative humidity (or dewpoint 
depression), wind speed, and wind direction were abstracted onto 


coding sheets from the above sources at all reported levels up to 


or 


abit B * o- var 
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and including the 200 mb level where possible. The abstracted data 
were then transferred onto punch cards for computer analysis. The 
height in geopotential meters Z of the desired levels was obtained 


using the standard equation: 


| 


pas ees 


*% « 
cn ;-RI In( Pax J (4.31) 
9.8 


in which i represents the ie level, R is the gas constant for dry 
a) iG sand lace is the mean virtual temperature in Kelvin degrees of the 
layer between pressures P. and epee At pressure levels where the 
temperature and moisture variables were not specifically given, they 
were obtained by linear interpolation between the reported values 
above and below that of the desired level. The interpolated values 
were used in (3.3.1). The calculated heights were used to obtain the 
wind speed and direction at the desired levels from the wind sounding, 
by linear interpolation where necessary. Missing data were not 
interpolated from surrounding stations. 

Several checks were made of the punched data to minimize 
the errors arising from the abstraction process. The punched data 
were checked manually against the coding sheets. The punched data 
were then listed by computer print-out and the print-out checked 
manually for apparent inconsistencies. Finally, the calculated heights 
of the 700, 600, 500, and 400 mb surfaces were checked against the 
heights reported on the adiabatic charts or RAOB computer listings. 
Height differences of six or more geopotential meters were considered 
sicomiticant forethisafinal check. 

Terrain heights were extracted at half-grid intervals over 
the geographical area bounded by grid points five to twelve inclusive 
in the y-direction, and three to eleven inclusive in the x-direction. 
The highest elevation within a 10 nautical mile radius of each point 
was used as the height at these points. To obtain the height at the 


full-grid points, the extracted terrain heights were averaged as 
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follows: for interior grid points 
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and for grid points on the boundary 


where Z ij is the terrain height at grid point (i,j). The 
Canadian terrain heights were extracted from World Aeronautical Charts 
ICAO, National Topographic Series (1:1,000,000). U.S. portions 

came from Topographic Maps of the U.S. (1:1,000,000) and National 
Topographic Maps (1:250,000), both published by the U.S. Department 

of Interior Geological Survey. The smoothed terrain contours are 


shown in Figure 7, and the Laplacian of the smoothed terrain is shown 


in Figure 8. 


3.4 The Objective Analysis 


The objective analysis program used was given by Glahn, 
Hollenbaugh and Lowry (1969). This was based on a program developed 
by Cressman (1959) which in turn was based on the general method 
described by Bergthérssen and Dodds (1955). It was chosen because it 
can analyze a height field without using the geostrophic wind 
relationship, although this can be included if desired (cf. Glahn and 
Hollenbaugh, 1969 ). Cressman (1959) points out that if wind is 


used in the analysis to obtain more accurate height gradients, then 
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Fig. 7 The smoothed terrain of the area for which the ageostrophic 
winds were analyzed. Contours are labelled in hundreds of meters above- 
sea-level. 
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Fig. 8 The Laplacian of the smoothed terrain of Fig. 7 using ae 
finite-difference interval of 185.3 km. Isopleth units are 1078 m-!, 


Ge a alee A a 7. 5 “a ; 7 
; a Blad? Sho. et, wal how we 
La Co wee - eed : : soipane : = : s 

~~ < ed 3% y”.4 Y aad 


i 
yiINgQoTIaCags eF2 io tae Velie aT" B84 So ciara’ et aehiee ete - + te Pa 


svodEe Av gar ; in. Bho briny ie bad lates, « =o. TUS. ee! One, cal ie se 


i : we : ; ce 


in areas of relatively dense observations the resulting analysis may 
differ appreciably from geostrophic if the wind and height fields do 
not agree well with the geostrophic approximation. The balance 
equation could be a suitable alternative since it involves the use 

of a stream function, and hence would try to force the resulting analysis 
to be non-divergent. However, inaccuracies in the observed data, and 
the necessity of solving a Poisson equation over the grid tend to 
produce an analysis which does not differ appreciably from that using 
the geostrophic relationship. To circumvent this problem and effect 
a considerable saving in computer time only the height field was 
analyzed and the geostrophic winds were obtained using equations 
C2543): 

The method of analysis consists basically of using the 
reported data to make successive corrections to an initial guess by 
making several passes (or scans) over the grid. In this study, the 
average value of the observed field was used as the first guess to 
initialize the grid. Bilinear or biquadratic interpolation is used 
to interpolate from the grid points to the station position. This 
occasions a change to the individual grid point values from all 
station observations within a prescribed radius of influence from the 
grid point. The corrections to the grid values can be made by three 


different methods: 
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in which Cc, 9 3 are the corrections to be applied, n is the number of 
Fae J 


Stations affecting the grid point, iF is the observed station 
d 


value, and AT is the interpolated station value. The weight function, 
b] 


W, is the one used by Cressman (1959) and is defined by: 


saul cual 
aioe 


where N is the radius of influence for the scan and d. is the distance 
th : ; ; : ; 
from the i Station to the gridpoint. The radius of influence can 


be varied for each scan. Corrections of type C, lead to a more rapid 


3 
convergence to the general level of the data within the influence 
rol ig el end We 

In regions of sparse data, there is a tendency for 
discontinuities to develop at the distance where the weight factor 


approaches zero (i.e. when d. approachesf{l). This necessitates the 


use of a smoothing operator. The one employed is given by: 


Neg ey 
oh = _b) i BAS 


and 


ee) 


, Fe See ia 2 os Ae 
b+t,y u-4, | b,jrt Cio! 


in which b is a number specified by the programmer and can be varied 
with each scan, 5. . is the smoothed value of the grid-point values 
> 


A, .. Clearly, if the smoothing parameter b is zero then no smoothing 
> 


will be done. 


It was necessary to use sixteen passes for the pressure- 
surface analyses and nine passes for the wind-field analyses in order 


to produce good agreement with subjectively analyzed samples. 
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3.5 The Analysis Procedure 

The raw data were analyzed as described in Section 3.2, 
with the pressure height and wind components being extracted at 25 mb 
intervals from 700 mb to 400 mb inclusive. A time-average of the 
heights for each level was used as the data input for the objective 
analysis to obtain the mean flow patterns. 

A geostrophic wind at each level for every station and time 
was obtained by objectively analyzing the height fields and using 
equations (2.1.3). This allowed Uy and vy to be computed. A time- 
average of the ageostrophic components was then calculated for each 
Station and level. The station values of UW, and V, were used for the 


A A 


Bellamy calculations. An objective analysis of Uy and V at each 


level was carried out to calculate the divergence and vorticity 


by the fundamental method. 
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CHAPTER IV 
THE RESULTS OF THE ANALYSES 


Hob biibntroduction 

In order to reduce boundary effects, and also because 
centered-space finite differences were employed, the ageostrophic wind 
components and the relative vorticity and divergence associated with 


them were computed for a reduced portion of the grid as shown in Figure 9. 
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Fig. 9 The geographical area covered by that portion of the grid 
for which the ageostrophic divergence and relative vorticity were calculated, 
Thin dotted line is the approximate position of the 1500 m ASL terrain contour. 


It should be noted that a reduction was not made over the 
eastern portions in order to be able to include the lee slopes of the 


Rocky Mountains. 
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The results of the divergence and relative vorticity 


calculations will be 


presented as vertical averages over the intervals 


700-625 mb, 600-525 mb, and 500-400 mb, which will be referred to as 


the low, middle, and 
purposes and to make 
overlays are located 
include the smoothed 


as shown in Figure/7 , 


high levels, respectively. Also, for comparison 
the charts more legible, several transparent 
inside the back cover of the thesis. These 
contours of the terrain for the reduced grid 


the average precipitation areas for the December 


and January periods studied and shown in Figures 5 and 6, and the 


terrain Laplacian of Figure 8. 


4.2 The Geostrophic Flow for the December Stud 


The time-averaged geostrophic flow patterns for three levels 
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Fig. 10 The time-averaged geostrophic flow of the December study 
for (a) 650 mb, (b) 550 mb, (c) 450 mb. Contour interval is 50 gpm. 
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The flow shows reasonably constant directions at all levels. 
Maximum geostrophic winds at all levels were found over the Washington 
area with values ranging from 20 m sah 
400 mb. 


at 700 mb to 30m ee at 
The average position of the upper low is indicated to be 


northwest of Vancouver Island. Confluent flow occurs along the coast 


to the east of Vancouver Island while diffluence is prevalent over 


the area east of 120°W. Weak ridging occurs at all levels in the 


vicinity of the Alberta portion of the Rocky Mountains. 


4.3 The Ageostrophic Wind Fields of the December Stud 


The direction of the time-averaged ageostrophic deviations 


showed little change with height as shown by the samples in Figure ll. 


© PRINCE GEORGE 


Sate oa a 


7 
\ 


Nes 


ee ek 
VERNON “ 
( 

i 


i290 [Om sec” 


Fig. 11 The time-averaged ageostrophic wind components of the 


December study for ‘a) 650 mb, (b) 550 mb, (c) 450 mb. 
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The strongest ageostrophic winds were found over Vancouver 
Island, reaching maximum speeds of 14 m sae at 400 mb. Comparison 
with Figure 10 shows that the ageostrophic component over Vancouver 
Island is predominantly across the contours toward lower contour values. 
Although the geostrophic winds show a general increase in speed at 
progressively higher levels, a similar trend is not apparent in the 


ageostrophic winds. 


4.4 The Ageostrophic Divergence Fields of the December Study 

The difference found by subtracting the divergence 
calculated by the fundamental method from that obtained by Bellamy's 
technique was obtained at the positions of the Bellamy triangle 
centroids for each level. The resulting 78 values had a mean of 
rae con sae with a standard deviation of Taeo,” ae A 
comment on these results is necessary. To compare the results in 
this manner is not strictly valid since the fundamental method and 
Bellamy method use different regions of data to compute the vorticity 
and divergence. In spite of this, a difference in sign occurred 
at only 14 of the 78 available points. If the Bellamy values were 
incorporated with the fundamental values by some method, the results 
would show only minor changes from those obtained by the fundamental 
method alone. Hence, only the fundamental results are presented in 
Figure 12 because of the more complete coverage obtained over the area 
of concern. 

The major convergence cell near 50°N, 125°W fits well with 
the previously noted confluent mean flow over this area. Similarly, 
the divergence area over Alberta is reflected by the diffluent nature 
of the mean flow over the eastern portions of the grid. The 
divergence cell centered near 46°N, 122°W is surprisingly strong 
considering the weak diffluence of the mean flow in this area. Gradient 
wind contributions can only account for about 20-25% of the observed 
values. Estimated gradient wind contributions over Montana and the 


-6 -1 
extreme southern portions of Alberta ranged from -2xlG_— sec in the 


low levels to ahi eau in the higher levels. The gradient wind 
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Fig. 12 The ageostrophic divergence fields of the December study as 
calculated by the fundamental method for (a) low levels, (b) middle 
levels, (c) high levels. Isopleths are in units of 107© sec- 
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cannot account for the strong nature of the secondary convergence cell 


which appears in the high levels over Idaho, although it is probably 


a contributing factor. The northward shift with height of the secondary 


divergence cell which is centered over 51°N, 115°W at low levels is 
accompanied by a weakening in the middle levels as it passes through 
the region of gradient wind induced convergence. A marked increase in 
the intensity of this cell is shown in the high levels as it moves out 


of the area of gradient wind influence. 


4.5 The Ageostrophic Relative Vorticity Fields for the December Study 


A comparison of the Bellamy and fundamental results by 


the method described in Section 4.4 gave a mean difference of 


“007107° FO a with a standard deviation of O7e107” ace As with 
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the divergence values, the difference of the vorticity values was 
primarily one of magnitude -- not of sign. Hence only the fundamental 


results are presented in Figure 13. 
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Fig. 13 The ageostrophic relative vorticity fields of the December 
study as calculated by the fundamental method for (a) low levels, 
(b) middle levels, (c) high levels. Isopleths are in units of 10-5 sec-!, 


When compared with the mean geostrophic flow of Figure 10, 
one finds that the ageostrophic relative vorticities have the opposite 
"pyolarity'' to that expected from the geostrophic motions. It should 
be noted that this difference would probably not be reflected in the 

4 


geostrophic vorticity field since it is of order yom sony: while 
the ageostrophic vorticities are of order lo tee 

A comparison of the low-level vorticity pattern with the 
terrain Laplacian of Figure 8 shows good agreement for the region of 
negative vorticity, but poor agreement for the area of positive 
vorticity. The middle and upper levels show little correlation with 
the terrain Laplacian. It is possible that the vorticity patterns, 
and by direct association also the ageostrophic deviations, arise in 
part from the mean geostrophic flow trying to conserve its absolute 


vorticity, since air moving northward would tend to gain negative 


relative vorticity. This would account for the stronger negative 
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vorticity values since the components of northward deviation are 


considerably stronger than the southward components. 


4.6 The Mean Geostrophic Flow for the January Study 

As indicated by Figure 14, there is little change in the 
mean flow patterns with height. 

At low levels there is confluent flow over Vancouver Island 
extending east to roughly 120°W, and generally diffluent flow over 
Alberta. The middle and high levels show a strengthening of the flow. 
The confluent zone of the low-level flow disappears at higher levels, 
but there is still some indication of diffluence over Montana and the 
eastern portions of Alberta. Maximum geostrophic wind speeds ranged 


from 25 m Sscuy at 700 mb to 45 m ae at 400 mb, 
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Fig. 14 The time-averaged geostrophic flow of the January study 
for (a) 650 mb, (b) 550 mb, (c) 450 mb. Contour interval is 50 gpm. 
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4.7 The Mean Ageostrophic Deviations for the January Study 


The mean ageostrophic deviations for the January case study 
for which samples are shown in Figure 15 do not exhibit the same degree 
of consistency with height found in the December case. 

The most notable change occurs in the vicinity of Vancouver 
Island. A basically geostrophic low-level flow is indicated by the 
relatively small westward ageostrophic winds in this area, A southward 
component appears over this area in the middle levels which increases 
in strength with height up to the 400 mb level. Other areas of the 
grid show the same basic features at all levels -- southward deviations 
over the eastern portions of Alberta, and east to northeastward 
deviations over Idaho and Montana. The center of small deviations near 
Edmonton in the low and middle levels shifts to the southwest at high 
levels to a position near Vernon, B.C. Small ageostrophic components 


are found at all levels over the southern portions of Washington State. 
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Fig. 15 The time-averaged ageostrophic wind components of the 
January study for (a) 650 mb, (b) 550 mb, (c) 450 mb. 
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The ageostrophic winds above Idaho, Montana, and the 
eastern portions of Alberta indicate a strong to moderate cross- 
contour flow at all levels. The strong southward ageostrophic component 
above 600 mb around Vancouver Island also gives strong crossecontour 
flow in this area, 

The cross-contour flow over Alberta and Vancouver Island 
is interesting in that it is directed towards higher geopotential, 


making it difficult to interpret as a frictional effect. 


4.8 The Ageostrophic Divergence Fields of the January Stud 


The fundamental results are presented in Figure 16. 


The ageostrophic divergence values of the Bellamy method show good 

- -1 
agreement with these results, having a mean difference of 1x10 sec 
with a standard deviation of 6x10°° iowa 


120 
a 


PORT HARD Bee cocks 


55 


CA ats 


>= 


“ “ 50 
GREfT i 
‘ ® Se 


10 a 


Bag, 168 


Fig. 16 The ageostrophic divergence fields of the January study as 
calculated by the fundamental method for (a) low levels, (b) middle 
levels, (c) high levels. JIsopleths are in units of 107 sec7l, 
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The low-level results show a center of convergence near 
50°N, 122°W and a center of divergence near 46°N, 118°W, with 
convergence prevailing over Alberta and British Columbia. The appearance 
of the southward ageostrophic component over Vancouver Island above 
600 mb results in a southwest shift of the main convergence center 
which reaches maximum intensity through the middle levels. The middle 
and high levels also show the appearance of secondary centers of 
convergence and divergence over Montana and Alberta, respectively. The 
weak curvature of the mean flow resulted in negligible gradient-wind 


CErectomatera le hevieslign 


4.9 The Ageostrophic Vorticity Fields of the January Stud 


The difference found by subtracting the fundamental results 


from the Bellamy results had a mean value of “i. @elGe sea with a 


standard deviation of Tic. The two sets of results showed good 


agreement in sign. The fundamental results are given in Figure 1/7. 
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Fig. 17 The ageostrophic relative vorticity fields for the January 
study as calculated by the fundamental method for (a) low levels, 


(b) middle levels, (c) high levels. Isopleths are in units of 107> sec7l, 
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The strong southward ageostrophic wind components in the 
middle and high levels over Vancouver Island are reflected in the 
northwest shift, and increase in strength at these levels of the positive 
center which appears in the low levels over the Central Washington area. 

While the low-level vorticity pattern shows some correlation 
with the terrain Laplacian of Figure 8, the agreement is poor in the 
upper levels. As with the December case, the relative vorticity of 
the ageostrophic deviations is of the opposite sign to that of the 


mean flow in the Same areas. 


4.10 Terrain and Condensational Heating Influences in the Results 

The main objective in undertaking this investigation was 
to try and detect terrain influence in the wind field at levels 
where direct interference by the terrain would not occur. Unfortunately, 
attempts to correlate the calculated ageostrophic deviations with the 
terrain features of Figure 7 were unsuccessful, with two possible 
exceptions. These exceptions involve the small area of precipitation 
located near 49°N, 122°W which occurred in both the periods investigated. 
An examination of the soundings for Port Hardy and Quillayute during 
the January period indicated latent instability in the levels near the 
surface which could be triggered by the forced ascent encountered by 
the air in these low levels as it surmounted the Coast Range. Air rising 
due to this instability could not penetrate much above 600 mb because of 
a strong inversion from this level to about 500 mb. The latent heat 
released by condensation in this rising air and the resulting precipition 
may account for part of the low-level convergence over this area. 
However, Aubert (1957) has shown that the largest values of convergence 
generated by condensational heating should occur below the level of 
greatest condensation. Hence, the observed maximum of convergence over 
this area in the 600 to 525 mb levels cannot be explained by either 
terrain or condensational heating effects. An analysis of the secondary 
center of precipitation for the December study produced similar results, 
with the exception that the factor inhibiting vertical motions was a 


lapse rate somewhat less than pseudo-adiabatic above 550 mb. 
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The soundings for the other stations in that portion of the 
grid for which ageostrophic winds were calculated showed that Port 
Hardy and Quillayute were the only stations for both periods with any 
significant degree of instability present. This implies that terrain 
effects were probably of a very minor nature at the levels analysed. 

The major precipitation area which occurred over Alberta 
during both the December and January periods was caused by weak frontal 
activity in the 850 to 700 mb layers and may account for the weak 
convergence in the lower layers over Alberta during the January period. 

It is apparent that some other controlling influence must 
be present that was not considered in Chapter II. The effects of 
horizontal thermal advection have been neglected to this point 


and will now be examined with regard to the results. 


4,11 Thermal Advection Effects for the December Study 

The December soundings for the stations enclosed by or 
boardering on the reduced grid of Figure 9 were examined for evidence 
of horizontal thermal advection. When the results of this investigation 
were compared with the computed ageostrophic deviations of Figure 11, 
it was found that in general, a southward ageostrophic component occurred 
with cold advection, while a northward component occurred with warm 
advection. In the light of these occurrences, the effects of a strong 
but slowly moving trough of warm air aloft which moved inland from the 
Pacific during the period are shown by the warm advection above and to 
the west of Vancouver Island. This warm advection is responsible for 
the strong ageostrophic convergence around Vancouver Island. The 
relatively light ageostrophic winds over Washington are a result of the 
time-averaging procedure. However, the strong divergence over this 
area reflects the dominant effect of cold advection following the passage 
of the trough of warm air. Warm advection in the middle and high levels 
accounts for the secondary cell of convergence found at these levels. 
The weak warm advection above 700 mb over Alberta and the Prince George 


area is off-set by strong advection of cold air below this level 


resulting in net ageostrophic divergence. 
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4.12 Thermal Advection Effects for the January Study 

As one might expect from the rather complicated pattern of 
the ageostrophic winds for the January study, thermal advection was 
somewhat more difficult to determine from the soundings. However, 
several features were apparent. The southward deviations over Alberta 
accompanied cold advection at all levels and resulted in the divergence 
which appears in the middle and high layers. As noted in Section 4.10, 
the low level convergence may be due to condensational heating effects. 

The soundings for Prince George and Port Hardy indicated 
that the strong southward ageostrophic component which develops over 
and to the west of Vancouver Island above 600 mb accompanies the 
advection of a warm dry air mass into this area at these levels, while 
weak cold advection was present below 650 mb. The narrow swath of 
ageostrophic winds having a northward component found over central 
British Columbia at all levels up to 500 mb was caused by strong warm 
advection below 500 mb in this area during the first half of the period. 
Thus, the convergence area which appears over most of British Columbia, 
Vancouver Island, and the northwest portions of Washington is a simple 


consequence of a rather complex thermal advection pattern. 


4.13 Summary of the Results 

The time-averaged ageostrophic wind components that were 
calculated show reasonable agreement with the characteristics outlined 
in a general study by Godson (1950). The divergence values calculated 
from the ageostrophic winds lie in the generally accepted range of 
107° Bet to Tome sec for synoptic scale atmospheric motions. The 
ageostrophic relative vorticities, with magnitudes of order oe sec. 
conform to the normal assumption that the divergent component of the 
wind field has a relative vorticity which is one order of magnitude 
less than that of the rotational (i.e. quasi-geostrophic) component 


- -1 
Which Le of order 10 a Secs 


There appears to be a weak contribution by condensational 
heating to the low level convergence found over Alberta during the 
period of the January study. Condensational heating may also contribute 


to the low level convergence found over British Columbia in both periods 
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Studied. Gradient winds had a small effect on divergence areas over 
Washington and Alberta during the December periods, but were negligible 
during the January period. 

It is felt that the apparent absence of terrain influence 
throughout the 700-400 mb layer during both periods may be caused by 
Ene stability of the pe mospiere over the area studied. The correlation 
between thermal advection and the regions of ageostrophic divergence 
and convergence suggests that horizontal thermal advection was the 


influence controlling the ageostrophic deviations. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS, WITH SOME 
SUGGESTIONS FOR FURTHER WORK 


5.1 Summary of the Methods Employed and the Calculated Results 


Two methods were outlined that are suitable for calculating 
the divergence and relative vorticity of the ageostrophic wind. The 
first method is referred to as the fundamental method and is based on 
the kinematic definitions of divergence and vorticity. The second 
technique is one developed by Bellamy (1949), and was intended primarily 
as a check on the values found using the fundamental methods. The 
merits and drawbacks of the two methods were discussed and it was 
decided to use the time-averaged ageostrophic wind fields for both 
methods in order to reduce the effects of observational errors. 

The ageostrophic component of the observed wind was 
obtained by subtracting an objectively analyzed geostrophic component 
from the observed wind. Possible sources of the ageostrophic wind 
component were discussed, including atmospheric compressibility, the 
gradient wind, diabatic heating by condensation, and terrain influences. 
Quantitative estimates of the contributions these factors could make 
to the ageostrophic divergence and relative vorticity were obtained 
where possible. A rough estimate was obtained of the uncertainty in 
the vorticity and divergence caused by errors in the observations, 
indicating that the values obtained by the Bellamy technique will 
generally be less uncertain than those of the fundamental method. 

Two case studies of the ageostrophic wind field were 
presented: one from December 1971, and the other from January 1972. 

For objective analysis purposes, and to facilitate the computations, 

a finite-differencing grid covering the geographical area studied was 
employed. The portion of the grid for which the ageostrophic divergence 
and relative vorticity were calculated included the southern portions of 


British Columbia and Alberta, the state of Washington, and portions 
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of Idaho and Montana. A smoothed topography for the terrain of these 
areas was obtained, and the Laplacian of this topographical surface 
was computed. 

Using the original radiosonde ascent recordings, the 
geopotential heights of the pressure surfaces from 700 to 400 mb at 25-mb 
intervals, as well as the corresponding wind vectors, were obtained at 
each observing station within the area covered by the grid, giving a 
total of thirteen levels. A time-average of the calculated heights for 
each level at each station was used to obtain the objectively analyzed 
mean flow for each of the two periods. The geostrophic wind was obtained 
for each station at every level and time using objective analysis and 
finite-difference techniques. This permitted the computation of the 
ageostrophic winds at the desired levels for each station and time. 

A time-average was taken of the ageostrophic components and used to 
calculate the Bellamy values of the ageostrophic divergence and 
vorticity. An objective analysis of these average ageostrophic 
components permitted the calculation of the fundamental values for the 
vorticity and divergence. 

The time-averaged ageostrophic wind fields for both periods 
conformed to the normal assumption that the divergent portion of the 
observed wind is generally an order of magnitude less than the rotational 
portion. Cross-contour components appeared at all levels over various 
portions of the grid. The sign of the Bellamy and fundamental results 
showed generally good agreement, with acceptable differences in 
magnitude. The ageostrophic divergence values occurred in the 


= -l 
: to 0 5 sec , while the 


re a) -1 
ageostrophic relative vorticities were of the order of -10 seq” 


synoptically acceptable range of 2107 


Condensational heating appeared to contribute to convergence 
in the 700 to 625—mb layer over Alberta during the January period. 
Diabatic heating may also have given some low level support to the 
convergence which occurred over British Columbia during both periods. 
Ageostrophic components arising from the gradient wind had a minor 
effect on the divergence fields during the December period, but were 


negligible for the January case. 
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No correlation was evident between the ageostrophic wind 
fields and the terrain, while the ageostrophic vorticities showed 
fair agreement with the terrain Laplacian in only the lowest levels 
studied. The ageostrophic winds and the divergence fields derived 
from them did show a good correlation with horizontal thermal advection 
that occurred during the periods studied, suggesting that this was the 


influence controlling the ageostrophic deviations for both periods. 


5.2 Conclusions 

No conclusions could be drawn on the nature of terrain 
influences in the upper flow during the periods studied. However, this 
does not mean that terrain effects were not present -- rather that 
the methods employed, the data available, and the prevailing atmospheric 
stabilities did not permit sufficient resolution of the wind field 
to be able to detect the results of terrain influence. The low level 
nature of aoe Gmail eee which occurred, and the calculated mean 
flow patterns resulted in only minor effects from the influences of 
condensational heating and the gradient wind. 

The rather unexpected ability of horizontal thermal advection 
to account for the dominant features of the ageostrophic divergence 
fields of both periods strongly suggests that this was the controlling 
influence on the ageostrophic wind fields. More work is required to 
clarify the relationship between thermal advection and the ageostrophic 
wind, However, it appears that the synoptic-scale ascent or subsidence 
accompanying thermal advection may be responsible for the non-geostrophic 


components of the flow for the two cases studied. 


5.3 Suggestions for Further Work 


Perhaps this investigation would be aided most by studying 
several more cases of strong flow across the Western Cordillera to 
try and determine if there is some feature of the ageostrophic wind 
field which appears regularly over some portions of the terrain. In 


this connection, and in the light of the results of the present study, 
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it would be advisable to check the atmospheric soundings closely for 
the prevailing stabilities and the levels at which condensational 
heating would be most significant. Also, it would probably be of 
some help in identifying terrain effects if the grid were expanded to 
include more area east of the Cordillera. This expansion would allow 
comparisons to be drawn between areas of strong and weak terrain 
effect. The grid expansion would probably prove most profitable over 
the Continental United States because of the data coverage available. 
It is felt that the use of time-averaged ageostrophic winds 
may have resulted in a loss of the detail necessary to identify terrain 
influences. Thus, the study of the data for specific times as opposed 
to time-averaged data may prove to be of some use in resolving terrain 
effects, although the results would be more susceptible to observational 
errors. It would also be interesting to use the data at specific times 
to examine the conditions in the ageostrophic wind fields during 
periods of lee cyclogenesis to see if there is some distinguishing 


feature which accompanies this phenomenon. 
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60 
LIST OF SYMBOLS 


Overbar (~~) denotes "time average" 


Asterisk (* ) denotes "actual value", excepting virtual 
temperature 


S denotes "uncertainty" 
Area as measured on a constant~-pressure surface 
The value of a variable at a grid point 


Station value of a variable interpolated from surrounding 
grid points 


Correction applied to a grid-point variable during an objective 
analysis 


Adiabatic speed of sound in air 

Specific heat at constant pressure of dry air 

Respectively, the divergence and vorticity on a reference surface 
Isobaric divergence and vorticity, respectively 


Isobaric divergence calculated by the fundamental and Bellamy 
methods, respectively 


Respectively, the isobaric divergence and vorticity of the 
gradient wind 


Distance on the grid from the a station 
= AOL gin 9 , the Coriolis parameter 
Gradient wind speed 

Acceleration of gravity 

Grid interval 


hist F 
Perpendicular distance from the i- vertex to the opposing 


side of a Bellamy triangle 

Curvature of a streamline or trajectory 

Length scale of horizontal velocity variations 
Length scale of vertical velocity variations 

Mass 

Radius of influence used in the objective analysis 
Frequency of oscillation 

Distance perpendicular to a streamline or trajectory 


The observed station value of a variable at position (x,y) on 
the grid 


The grid positions of an observing station in the x and y- 
directions, respectively 


The grid positions of the principal grid point 
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Pressure 


Isobaric relative vorticity calculated by the fundamental and 
Bellamy methods, respectively 


The specific gas constant for dry air 
Radius of curvature of a streamline or trajectory 
Radius of the earth 


The image-plane distance from the North Pole to a point on the 
grid measured along a meridian 


Specific entropy 
Distance along a streamline or trajectory 
The grid-point value of a variable after a smoothing operation 
Temperature (°K) 

( apc) T, the virtual temperature (°K) 

l+w 

Time 
Horizontal scale velocity 
Observed wind components on an isobaric surface 
Calculated ageostrophic wind components on an isobaric surface 
Calculated geostrophic wind components on an isobaric surface 
Observed wind vector 
Ageostrophic wind vector on an isobaric surface 
Geostrophic wind vector on an isobaric surface 


N'-d, , a weight function 


Nerds 
i 


Scale vertical velocity 

Mixing ratio 

Terrain height 

Geopotential height of a pressure surface 

The direction of the ageostrophic wind component 


Angle between the reference meridian and the meridian through 
a station 


Angle between the reference meridian and the meridian through 
the principal-grid-point 


th 
Wind direction at the i- vertex of a Bellamy triangle 


Wind direction measured clockwise from some reference direction 
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Lsobaric coefficient of thermal expansion 


atl 2 
Azimuth of the side opposing the i vertex of a Bellamy 
triangle 


Finite difference 
Contour spacing 
A measure of the temperature difference in a fluid 


oo , thermal diffusivity of dry air 


c 
P 


Dynamie viscosity 

Kinematic viscosity 

stability factor 

Density 

Variance 

Latitude 

Angle between the perpendicular to a contour and the x-axis 
Co-latitude 


= Fo ,» velocity perpendicular to an isobaric surface 
dt 


-t3 a 57 kpg3e , the three-dimensional del-operator 


The two dimensional del-operator on a constant-pressure 
surface 


Two-dimensional Laplacian on a constant-pressure surface 
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APPENDIX [I 


DIVERGENCE AND VORTICITY USING THE BELLAMY TRIANGLE METHOD 


Consider the equation of continuity in the form: 


ens (hada ts 


— 


a 
0 
<4 


al 
cg 


as applied to a parcel of air moving on a constant pressure surface 
having an area on this Surface of A, unit thickness, and containing 


Mass m of air, Applying (A.1.1) to this parcel would yield: 


-A ad (p)=ve¥ 


or 


V (BD) 


dA = V 9 
Ata bee 


Bais 


UNE aah 


or 
Dee ede Gate, 
A dt 
where D_ represents the velocity divergence of the parcel as it moves 
Pp 
Let the surface area of the parcel be 


along the pressure surface. 
represented by a triangular shape as shown in Figure Al. 
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Fig. A.1 The geometry of the Bellamy Triangle. 


The area of the triangle is given by 
Vetoe ea tome 
és 


If we assume that the wind speed is zero at vertices b and c, 
and increases linearly to ve at vertex a, then the outflow (or inflow) 
of air from the volume of unit thickness in some arbitrary unit time 


At is represented by 


AA=1(h+n.) be -1 h, be 


whence 


(A.1.4) 


e 


aoa: 
At 


ha 


p 
Ba 


where D represents the partial divergence at vertex a. Referring to 
a 


Figure Al, one finds that 


7, = —V,$in Gye 7) 
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Thus (A.1.4) can be written: 
eae ONIN Ce 
ha 


Carrying through a similar argument for the remaining 


vertices and summing the results gives: 


IDEs = 32 Me sin (a a) Cais) 


yep L 


where the numbers 1,2,3 are substituted for the letters a,b,c and BY 
represents the azimuth of the side opposing vertex i. Results of this 


operation yield positive values for divergence and negative values for 


convergence. 


To obtain the partial vertical component of vorticity at the 


i vertex one need only add 90° to the wind direction and use (A.1.5). 
Zhus: 


and 


3 
a Ve-eus Fs - Ble 
Q d, a (Ai- &; ) (A.1.6) 


Use of (A.1.6) yields positive values for cyclonic vorticity 


and negative values for anticyclonic vorticity. 
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APPENDIX ITI 
GRADIENT WIND VORTICITIES AND DIVERGENCES 


Employing the prime notation to indicate differentiation 


with respect to x, the radius of curvature of a streamline (R) is 


given by: 


y, 
ee escyS) a 


CAT 24) 
4” 


Since G is assumed constant, equations (2.3.4) become: 


CAL 2) 
2 an 


£4 
ry) 
v 


C2) 


also, since oe pO) > we Lind: 


38 = dp dx (A.2.4) 
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bz=arctan (y’) (A.2.5) 
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Fig. A.2 The finite-difference elements for evaluating the partial 
derivatives in the equations for the gradient wind divergence and relative 
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Referring to Figure A.2 and using finite differences to represent the 


partial derivatives: 
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and using (A.2.7) gives: 


= = SIN 8 


Using (A.2.8) and 


eens 
an R 


Thus; (A.2.2) and 
/ 
ee 
R 


Qye 


where R is given by (A.2.1). 


Livy?’ 


(A.2:8) 
(A.2.6) we obtain: 
ae Sin B 
(A.2.3) become: 
(A.2.9) 
(A.2.10) 


Equations (A.2.9) and (A.2.10) were evaluated on an 


-1 
IBM 360/67 computer using G=40 m sec , and x and y were taken to be 


measured in grid intervals of 


185.3 km. Values were calculated for 


x=-10 to 10 in increments aps MO ae 
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Since the parameters (i.e. R, y', etc.) do not change if a 


ditterent. value ot G is chosen; i1£ is clear that: 
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where Dar and Qa are the divergence and vorticity estimates for 


the new gradient wind speed G'. 
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APPENDIX [LE 


DERIVATION OF THE UNCERTAINTIES IN THE DIVERGENCE AND 
VORTICITIES CALCULATED BY THE FUNDAMENTAL AND BELLAMY METHODS 


In finite difference notation, the geostrophic wind components 


calculated from the height field are given by: 


(A731) 


Mey aH Vas fie 


where H is the grid interval, the image factor of the map has been 


neglected, and the subscripts refer to the points as given in Figure A.3. 


e 
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i) @ 0 ® 3 
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Fie, A.3 The grid-point ordering for calculating the geostrophic winds 


on a constant pressure surface. 
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In the following derivations, only the equations for the 
west-east component of velocity (u) will be presented where possible, 
since those for the north-south component (v) will have the same forn. 

In that the pressure heights at the grid points are just 
extrapolations from the data points, the uncertainties in the initial 
data set (i.e. the heights calculated from the radiosonde data) will 
be communicated to the grid point values. To try and calculate how 
much of the initial uncertainty in the height values will be transferred 
to the grid points by the objective analysis routine would be difficult. 
In lieu of this, it was assumed that the uncertainty ( SZ Vai et he 
initial data set would apply equally to the grid values. 


Using this assumption, (A.3.1) can be expressed as: 


ike = wg © g 5H Gis) 


where u_, v_ are the calculated values, and superscript * is used to 
8 & 
indicata “actual value"; 


Using (A.3.2), the ageostrophic components of the wind are: 


ux = (Gi By) WE OS eg ae) (rae3) 


Se a) 


where u and v are the observed wind components. . Letting Naame 


(ATS. 3) can, De written as: 


a= Ut (Surg $2 ) (A.3.4) 
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The average of the ageostrophic components at a given 


level at a given point can thus be expressed as: 


AT Gn Uys I : ( Sut Sane (At 3-5) 
A a dharare 


where L is the number of values over which the average was taken. 


However, since SL and SZ are not considered functions of time, 


(A.3.5) can be simplified to: 


5% asia. ieee 


a BSG el yi ne 
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Employing finite differences to evaluate the derivatives 


gives the average fundamental divergence ( ies Ne 


ie ~ Dea Oise) aq as 4 Cee 


The uncertainty involved in the fundamental vorticity 


would be identical. Thus: 


SN Oh =| ahs Cael 2g 2% | (A538) 
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where ae and $Q, are the uncertainties in the average divergence 
and vorticity calculated by the fundamental method. 
Let us now examine the Bellamy method of divergence 


calewl ation 4s. e1ven. oy (1.2.2) 


D.= “¥ sin (g-«) (A.3.9) 


Since the Bellamy calculations utilized the average 


ageostrophic components, the wind speed in (A.3.9) is given by: 


V= [ie ss ee CAB BELO) 


This means that the angle c{ represents the direction of 
the average ageostrophic wind as measured in a clockwise manner from 


the y-axis of the grid, or: 


oo = arcton| o } (A.3.11) 


From calculus: 


AINE sets, a Saye 
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mek 


and using (A.3.10) gives: 


$V= | (tie sateeee i Say ) (A.3.12) 
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Similarly 


Mi = ots OLE See Sar 
du ov 


and using (A.3.11) gives: 
(CARS 13) 
Now, from (A.3.9) 


oD,=- sin ae oV ty cos (p-«) Sa (A,3.14) 


Substituting for V, $V, and Sa and simplifying gives: 


SDE 2 (cos (p- “)[%,Sa,- ASU] 


| 
h yh + Ty 
~ sin(p-«) pie. Eee sv} 


However, since (A.3.9) must be evaluated at the three 


vertices of the triangle and the results summed, the total uncertainty 
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in the Bellamy divergence becomes: 


(AGS 25) 


Carrying through a similar analysis for the Bellamy 


vVorticities gives: 


SGee = 2 i) io [G, sty 


hla + Pe 


pS 2 00 ay oe ee (A.3.16) 
+U,Sv,|+ sin(a-«)[¥, Su, 


v.13]. 


Making the following approximations: Hu 3 oe Sue Sy, 


and taking $1 (6-x) , COS (g-«) to be of order one will simplify 


Cameo ES) cand (A,3.16) "tot 


SD 
oh ae 


OSes 22 é | Gee 


ae ee! comag erst 


[kee 8 


a) 


+ 


— | S 
ae BY le G09} + 
ee ee as cae 


<’.@ , 
,i6 (db! & Sng ~| Ege 
F if r * 
; -. . * 
; 
i* 
- é | j iyi Cpe? & te 
per! s tai2lisisywVv 
. & 7 A. 
ut | ; ) i = ee 
.¥ 
at , v - 
yy (ea UC ee & e és 4 
7 . =| 
A 
‘WE ~we 7 ’ lie? eit? gna 
eUr rapes pa! he § nt (aglaes eg 1700 gnige? Soe 
ps 000,49 tee (ELE 
(55,4, A> k e i¢ ae : 7 faa G ig 
| =f A\ AG ) pe = 
M4 Ay ° ag? 
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where it has been further assumed that the terms involving itv and u, 
have some average value for the three vertices. 


With the above assumptions (A.3.8) becomes: 


SD e $Q.+ ia (du+ g 2 | (A.3.19) 


Using CAl3517)> (A¢3218) and -(A.3519) “tosevaluate the 
ratios oa Sie and 2, y, S Q; gives the result 
that the uncertainties in the Bellamy divergence will be less than that 
of the fundamental method if aE @) while thessame result 
applies for the vorticities if a> Q 
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APPENDIX IV 


THE EQUATIONS FOR TRANSFORMING GEOGRAPHICAL 
POSITIONS INTO GRID POSITIONS 


Let the grid point (15,1) be called the principal-erid-peint., 
Let the grid position of the principal-grid-point be P 
Let the grid position of the station be Eee 


Then the grid positions for the station will be given by: 
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x =80-6 
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6, 80° W 


Fig. A.4 The image-plane geometry for transforming latitude and 
longitude into grid-point positions. 
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Referring to Figure A.4 one finds that: 


Poy ty (7 Orn) ~ 7 Ce) cos 5 | 
(A.4.2) 


Kae {rye Sin “en 


in which r is the distance in kilometers from the pole to the point as 
measured along a meridian on the image plane, Se is the reference 
co-latitude (40°N), The is the co-latitude of the principle-grid-point 


(40°N). 


simi barly: 
ag jes) AC OIS = 
(A.4.3) 


x 


i i { (Ys) sin %s | 


where We is the co-latitude of the station. 


For a polar stereographic projection, Saucier (1955) gives 


the expression for r(‘ ) as: 


PIAS I Ce COE jah (A.4.4) 
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in which Y is the co-latitude of the desired point, Ro is the 
radius of the earth (6371.22 km), and Ve is the standard co-latitude 
for the map projection (30°N). Equation (A.4.4) incorporates the image 


as 
scale meltcos Yo . 
ltcos ¥Y 


Combining (A.4.1), (A.4.2), (A.4.3) and using (A.4.4) 


gives: 
yee lo + FT ( tan Ye 6O$ a» -tan Ys cos os ) 
° 4 2. 2 
(AR 4.5) 
pene ele (tan 7s sina, - tan Ye sin %p) 
H 2 x 
where 
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TERRAIN LAPLACIAN 


DECEMBER PRECIPITATION AREAS 


JANUARY PRECIPITATION AREAS 


